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ABSTRACT 

Propagation of radio waves in the ultrarelativistic magnetized electron- 
positron plasma of pulsar magnetosphere is considered. Polarization 
state of the original natural waves is found to vary markedly on 
account of the wave mode coupling and cyclotron absorption. The 
change is most pronounced when the regions of mode coupling and 
cyclotron resonance approximately coincide. In cases when the wave 
mode coupling occurs above and below the resonance region, the 
resultant polarization appears essentially distinct. The main result 
of the paper is that in the former case the polarization modes be¬ 
come non-orthogonal. The analytical treatment of the equations of 
polarization transfer is accompanied by the numerical calculations. 
The observational consequences of polarization evolution in pulsar 
plasma are discussed as well. 
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1 INTRODUCTION 

1.1 Empirical model of pulsar polarization 


The radio emission observed from pulsars is typically characterized b y a high percentage of linear polar- 
ization. Within the framework of a well-known rotating vector model liRadhakrishnan fc Cookelll969 ). the 
orientation of pulsar polarization reflects the magnetic held geometry in the emission region and there¬ 
fore the position angle (PA) changes monotonically as the pulsar beam rotates with respect to an ob¬ 
server. The characteristic S-shaped swing of PA across the pulse is inde ed observ ed in a number o f pul¬ 
sars. In addition, PA may show abrupt jumps by approximately 90° (e.g. iManchester. Tavlor fc Hugueninl 
Il97.t . testifying to the presence of the two orthogonally polarized modes (OPMs). The early studies 
of this phenomenon have reve ale d that for each of the OPMs PA roughly follows the predictions of 
the rotating vect or model feac ker. Rankin fc Ca.rrmhelllll 97ft and mode changing is a stochastic process 

Further on t he OPMs have been rec ognized as a fundamental feature of 

A comprehensive analysis of the 


the rotating vector model (Backer, Rankm & Campbell 1976) and mode 
JCordes. Rankin fe Backeilll978t). Further on the OPMs have been recog n 
pulsar radio emission IIBacker fc Rarikirilll 98(1 Istinebring et al . 1 984aflt . 
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2 S. A. Petrova 


observational data has proved an idea of superposed OPMs: At any pulse longitude the radio emission is 
believ ed to present a n i ncoherent mixture o f the two OPMs, whose intensities vary randomly from pulse to 
pulse iMcKinnon fe Stinebrinelll998l. l 200cl) . 

The plasma of pulsar magnetosphere does allow two types of non-damping natural waves, the ordinary 
and extraordinary ones. They propagate in a superstrong magnetic field, generally at not a small angle 
to the field lines, and, correspondingly, are linearly polarized in orthogonal directions. The electric vector 
of the ordinary wave lies in the same plane as the wavevector and the ambient magnetic Held, while 
the extraordinary wave is polarized perpendicularly to this plane. The origin of th e two types of natural 
waves is attributed either to the two distinct radio emission mechamsnr^JagjMi •Kiiiiion Tour : or to the 
partial conversion of the ordinary mode i nto the extraord inar y on e iPetrovalbOOIf) . A recently discovered 
anticorrelation of the OPM intensities (Ipdwa.rds fc StannerdliiofiiFl favours the latter scenario. 

It should be noted, however, that direct identification of the observed superposed OPMs with the 
natural modes of pulsar plasma faces serious difficulties. First of all, some circular polarization is always 
present in pulsar radiation. Usually it is much lower than its linear counte rpart but not negli gible. It has 
been noticed that the sign of circular polarization is well correlated with PA llCordes et al.ll978h . that is the 
two OPMs have circular polarization of opposite signs and can be regarded as purely orthogonal elliptical 
modes. As for the theoretical interpretation, immediate switching to the case of elliptical natural waves 
jMglros^^^tc rnehamlll97^: iMelroseJ ll97g; Ivon Hoensbroech. Lesch fe Kunzll'1998; iMelrose fc Lu3l2004at: 
ILuo fc Melrosell2004l) seems problematic. The ellipticity of the natural waves may result from the gyrotropy 
of pulsar plasma (caused by difference in the distributions of electrons and positrons), if only the waves 
propagate quasi-longitudinally with respect to the magnetic field. Although the plasma gyrotropy is very 
probable, the regime of quasi-longitudinal propagation is not characteristic of pulsar magnetosphere. It can 
be the case only at some specific locations (most likely, close to the magnetic axis, where the divergence 
of magnetic field lines is less significant) and cannot account for the elliptically polarized natural waves 
observed throughout the pulse and over a wide frequency range. 

Recent thorough studies of the single-pulse data have introduced further complications into the picture 
of pulsar polarization. It has been found that the observed fluctuatio ns of the Stokes parameters cannot be 
expla ined solely by the pulse-to-pulse variation of the OPM intensities iMcKinnonkOQlrlEd wards fc Stanoerj 
120041 ). To account for the observations it has been sugge sted to complement the OPMs with a randomly po¬ 
larized component of unknown nature JMcKinnonll2004ll . Alternatively, the same results can be interpreted 
as a consequence of pulse-to-pul se jitter of both the ellipticity and PA dKarastereiou. Johnston fc Kranieil 
l200.dl : lEdwards fc Stanner jl2004lf . However, even this generalized picture based on the OPMs with the ran¬ 
domly varying v ector of the Stokes parameters appears incomplete. In some cases, th e modes are clearly 
non-orthogonal I Karastergfon et al .11200(1 iMcKinrionlhoftil : iF.d wards fc Stannerdtefiiull . It should be men¬ 
tioned that the non-orthogonality of polarization modes manifests itself not only in PA, but also in the 
circular polarization. In particular, the same PA may be accompanied by the circul ar polarizatio n o f any 
sense , the correlation between PA and V being less perfect at higher frequencies llKarastergiou et al.ll200ll 
(2003). 

Diverse and complicated behaviour of the single-pulse polarization as well as its strong frequency 
dependence motivate our study of the propagation effects in pulsar magnetosphere. 


1.2 Wave mode coupling in pulsar magnetosphere 

Pulsar radio emission is believed to be generated deep inside the tube of open magnetic lines. Then it 
should propagate through the flow of an ultrarelativistic highly magnetized electron-positron plasma. As a 
result of propagation effects, polarization of the radio waves may evolve significantly. In the vicinity of the 
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Polarization transfer in pulsar magnetosphere 3 


emission region, the ordinary and extraordinary waves are linearly polarized in orthogonal directions and 
the plasma number density is so large that the geometrical optics regime holds: the natural waves propagate 
independently, with the electric vectors being adjusted to the orientation of the ambient magnetic field. As 
the plasma number density decreases with distance from the neutron star, the difference in the refractive 
indices of the waves decreases as well, and finally the scale length for beating between the modes becomes 
comparable to the scale length for change in the plasma parameters. Then the polarization planes of the 
waves have no time to follow the local magnetic field direction, geometrical optics approximation is broken 
and wave mode coupling starts. Typically this occurs in the outer magnetosphere, at distances of a few 
tenth of the light cylinder radius. In the region of wave mode coupling, each of the incident natural waves 
becomes a coherent sum of both natural waves peculiar to the ambient plasma, with the amplitude ratio 
and phase difference varying along the trajectory. Correspondingly, the ellipticity of the waves increases 
with distance and the major axis of polarization ellipse is monotonically shifted, so that it no longer reflects 
the orientation of the ambient magnetic field. Further on, as the plasma density decreases considerably, the 
waves decouple from the plasma and propagate just as in vacuum, preserving their elliptical polarization. 
Therefore the process of wave mode coupling is usually called polarization-limiting effect. 


This effect has long been used to explain the origin of c i rcular polarization in pulsa r radio emis¬ 
sion lichens fc Rudermanlll97fll : iR.adhakrishnan fc RankiiilllQflrt : iLvubarskii fc Pet,rovalll999ll . The numer¬ 
ical tracings of polarization evolution in pulsar plasma have d emonstrated that the mode coupling effect 
is stro ng enough to have marked observational consequences I Petrova fc LvnbarskiilEoont IPetrovalboinl 
l2003alf . 


Polarization evolution of radio waves in pulsar plasma differs significantly from the evolution in the 
interstellar medium. In contrast to the case of Faraday rotation, within the pulsar magnetosphere the natural 
waves propagate quasi-transversely with respect to the magnetic held and have linear polarization. This 
rather correspo nds to the Cotton -M outon birefringence (or so called generalized Faraday rotation, in terms 
of the paper bv lKennet fe Melrose! Il998 j. Another important distinctive feature of polarization evolution 
in pulsar magnetosphere is that it takes place in an essentially inhomogeneous medium. The magnetic held 
of a pulsar has approximately dipolar structure, and furthermore, because of continuity of the plasma how 
in the tube of open magnetic lines, the plasma number density decreases rapidly with distance. Thus, the 
character of polarization evolution in pulsar magnetosphere is quite specihc, though the underlying physics 
of birefringence is certainly the same. 


As a result of the mode coupling effect, the outgoing waves acquire purely orthogonal elliptical polar¬ 
ization, matching the empirical representation of superposed elliptical OPMs. It is important to note that 
the degree of circular polarization of the modes and their shift in PA are related to each ot her, bot h being 
determined by the parameters of the plasma flow in the region of wave mode coupling dPetroval!l2003aj. 
Hence, the observed pulse-to-pulse variations in the ellipticity and PA of the OPMs can be attributed to 
the fluctuations in pulsar plasma. Besides that, the propagation origin of pulsar polarization should imply 
a correlation between the values of the ellipticity and PA of the OPMs at a given pulse lon gitu de. An_evi- 
dence for such a correlation has recently been found in iFdwardd ll2004h (for more details see IPetrovalboQfiil . 
Thus, the mode coupling effect can account for a number of important features of the observed single-pulse 
polarization. 


At the same time, the question as to the origin of non-orthogonality of polarization modes still remains 
open. The o bservational manifestations of this phenomenon have recently been repo r ted in a number of 
papers (e.g . lEdwards fc Stappersl I 20 q 3 : iMcKinnonl I2OO4I : iR.amachandran et alJ I2OO4I : iKarastergiou et alJ 
12003 . 1200In . In the present paper, we concentrate on a more detailed treatment of polarization evolution in 
pulsar magnetosphere, which, in particular, explains non-orthogonality of the modes. 
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1.3 Statement of the problem 


Polarization evolution in pulsar magnetosphere has previously been considered in the approximation of a 
superstrong magnetic held. It means that in the rest frame of the plasma flow the radio frequency, u >', is 
much less than the electron gyrofrequency, u>h- In other words, the radius of cyclotron resonance, where 
u / = u>h, has been assumed to be infinitely large. At the conditions relevant to pulsar magnetosphere, the 
radius of cyclotron resonance is often somewhat larg er t han t hat of the mode coupling region, but generally 
these quantities are of the same order of magnitude fearnardlll98fil . see also equation (15) below). Therefore 
it is reasonable to inspect the role of the cyclotron resonance in the evolution of pulsar polarization. 


In application to pulsars, the cyclotron absorption h as been considered in a number of papers llBla.ndford fc Scha.rlemanrl 
ll97dlLvubarskii fc PetrovJl994IPetrovJ200ll2003hll and found efficient, especially in case of small pitch- 
angles of the absorbing particles. In these papers, it is assumed that the resonant photons interact with a 
system of absorbing particles rather than with the plasma. The main motivation for such an assumption 
is that the resonance region lies in the outer magnetosphere, where the plasma number density is small 
enough. Indeed, in case of pulsars, taking into account the plasma effect on the process of cyclotron ab¬ 
sorption introduces only small corrections to the absorption coe fficients of the natural modes and does not 
change the total intensities of outgoing radiation considerably l lLvubarskii fc Petrovalfl998h . At the same 
time, cyclotron absorption in the plasma may markedly affect polarization evolution of radio waves. The 
contribution of cyclotron absorption, though not large quantitatively, may appear comparable to that of 
the mode coupling effect, modifying the final polarization of a pulsar drastically. 


Given that the regime of geometrical optics is still valid in the region of cyclotron resonance, the 
ordinary and extraordinary waves are absorbed independently, with the absorption coefficients being slightly 
different. Note that this difference is purely the plasma effect, not characteristic of a simple system of 
absorbing particles. Beyond the resonance region, the waves propagate in the weakly magnetized plasma 
with rapidly decreasing number density and finally suffer the mode coupling. It should be noted that in case 
of weakly magnetized plasma the limiting polarization differs substantially from that for the superstrong 
magnetic held. 

Given that the natural waves pass through the coupling region before the cyclotron resonance, in the 
resonance region each of the waves presents a coherent mixture of the two natural waves. Since for these 
constituents absorption is not identical, the wave polarization changes considerably. It is important to 
note that for the two incident waves polarization evolution is not the same, since they contain different 
portions of the ordinary and extraordinary waves. As a result, polarization states of the outgoing waves are 
non-orthogonal. 


In the present paper, we concentrate on the analytical consideration of polarization transfer in pulsar 
plasma, taking into account the effect of cyclotron resonance. The plan of the paper is as follows. In Sect. 
2, the main equations are derived, which describe the evolution of the wave fields in the inhomogeneous 
hot plasma embedded in the magnetic field. The basic numerical estimates are also given there. In Sect. 3, 
we solve the equations of polarization transfer in the two limiting cases, when the resonance region is well 
below and well above the mode coupling region, respectively. Section 4 contains the results of numerical 
tracing of polarization evolution. In Sect. 5, the observational consequences of polarization transfer in 
pulsar magnetosphere are discussed. Section 6 contains a brief summary. A statistical model of singl e-pulse 
polar ization based on the propagation effects studied will be developed in the forthcoming paper llPetroval 
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2 GENERAL THEORY OF POLARIZATION EVOLUTION IN PULSAR 
MAGNETOSPHERE 

2.1 Basic equations 

Let radio waves propagate in the ultrarelativistic highly magnetized electron-positron plasma of pulsar 
magnetosphere. The probl em on polarization evo l ution of t he waves in case of infinitely strong magnetic 
field has been considered in lLvubarskii V Petroval ill a.nd lPetrova V Lvnharsklil 1200(11 . The wave prop¬ 
agation in a cold relativistically streaming plasma embedded in the magnetic field of a finite strength has 
been studied in IPetroval ( 2001 1. However, in that paper the plasma number density has been assumed low 
enough for the polarization evolution to cease well below the radius of cyclotron resonance and the wave 
passage through the resonance has not been considered. Below we derive the generalized equations which 
allow for the cyclotron resonance and are applicable to the more realistic case of a hot plasma. 

The evolution of the wave fields E and B is described by the Maxwell’s equations: 


ILO 

V x E = —B. 

c 


_ 4-7T \ - . 

V x B =-E -|— y j a . 

c. c. Z—J 


(1) 

( 2 ) 


Here j a is the linearized current density of the electrons or positrons, summation is over the particle species, 
and the time dependence is taken in the form e -1 " 4 . We are going to trace polarization evolution starting 
from the place where refraction is already inefficient and the natural waves propagate straight. Hence, one 
can choose a three-dimensional Cartesian system with the z-axis along the wave trajectory and the xz-plane 
being the plane of the field lines of the dipolar magnetic field of a pulsar at the starting point. Then all the 
quantities in equations (l)-(2) depend only on the z-coordinate. 

Since the refractive indices of both natural waves, the ordinary and extraordinary ones, are close to 
unity, the wave electric field can be presented as 


Ex.V.Z - Ex 


where E x , y ,z vary weakly over the wavelength: 


d E x 


Az 


E — 

J-Jx.y.z 


(3) 


(4) 


All the other quantities entering equations (l)-(2) can be presented similarly. Further on the tildes will be 
omitted. The representation (3)-(4) corresponds to the case of a weakly inhomogeneous medium. Typically 
the scale length for change in the plasma parameters is about the altitude above the neutron star, 10 7 — 10 8 
cm, and the assumption of weak inhomogeneity is justified. Note that equation (4) requires also the length 
of the cyclotron resonance region to be much more than the wavelength. This is valid for the hot plasma, 
in which case the wave passes successively through the resonance with the particles of different energies. 
Substituting equation (1) into equation (2) and making use of equations (3) and (4) yields: 

AE X , 27 t \ ' . 

TTT + — 


d E v 27r v—i 

H7 + t5>»= O’ 


( 5 ) 


^ 47 ri . 

E z -\ -N jaz — 0 . 

UJ ' 


© 0000 RAS, MNRAS 000, 000-000 









6 S. A. Petrova 


As can be seen from the above equations, E z -C E X:V , that is the waves are practically transverse; hereafter 
it is taken that E z = 0. Furthermore, for the components of jo, in the first two equations it is sufficient 
to use the expressions obtained in the approximation of a homogeneous medium, where the fields and 
currents are oc e b ~ z , since the corrections for the medium inhomogeneity are too sm all (cf . e quati on (4)) . 
The conductivity tensor of a hot magnetized homogeneous plasma is well known (e.g. lMikhailovskvin~97fill . 
however, it is typically written for the coordinate system with the z-axis along the magnetic field direction 
rather than along the wavevector, as is necessary for our problem. The corresponding transformation of the 
conductivity tensor is performed in Appendix. Using those results in equation (5), we find finally: 

^ ^ [Ab x {E x b x 

a z 2 c 

d Ey IUJ r , _ , 

~dT + 2 c \- Ab v( E * b * 

Here b x ,y,z are the direction cosines of the ambient magnetic field, 


+ Eyby ) — BE X + iGEy] — 0 , 
+ Eyby) — BEy — iGE x ] = 0 . 


( 6 ) 


^E 


*/<* (7f>) 


2 

UJ H 


-y a iOa 


B = 


Wpa/ a (7a)7a (1 ~ P a b z ) 2 


G = 




^pa fa('y<*)(q a /e)(u H /u) [Poe - b z 


(7) 


— / 4-7re 2 r; 


is the plasma frequency, n the number density, u/ = 017(1 — Pb z ) is the wave frequency in the 

frame of the plasma particle moving at a speed of p = v/c, 7 = (l — /3 2 ) 1//2 the Lorentz-factor, q a = ±e, 
/(q) the distribution function of electrons or positrons with the normalization f f{ q)dq = 1 , the double 
square brackets stand for the Landau integral: 


F( 7 ) 


: v.p. 


F(7)d7 


+ ni 


5(ujh — o ; , ) l ? ( 7 ) d7 

UJH + oj' 


( 8 ) 


the first of the above integrals is taken in the principal value sense. 

Equations ( 6 )-( 8 ) describe polarization evolution of the waves in the hot magnetized weakly inhomo- 
gene ous p lasma allowing for the effect of cyclotron resonance. These equations coincide with equation (12) 
in IPetroval (l200lh in case of cold plasma, with the distribution function /(q) = <5 (7 — 70 ), if only a part of 
the wave trajectory well below the resonance radius is considered. 

It should be noted that equation ( 6 ) incorporates the current density in the limit of small pitch-angles, 
ip, of the plasma particles (for more detail see Appendix): 


u' ip 
loh 6 


< 1 , 


(9) 


where 8 ss b 2 + by is the angle between the wavevector and the ambient magnetic field. In the resonance 
region, equation (9) means that ip 6. Since in pulsar magnetosphere the optical depth to cyclotron 
absorption is typically large (cf. equation ( 12 ) below), the abs orbing particles can increase their p itch- 
angles significantly and the approximation ip 8 can be broken llLvubarskii fc Pet roviJlTinTj : IPet.roval2002i 
l2003l Jl . This especially concerns low enough radio frequencies which meet the resonance condition at higher 
altitudes above the neutron star and thus interact with the particles whose distribution function has already 
evolved because of absorption of the higher-frequency radio photons. At the same time, absorption by the 
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Polarization transfer in pulsar magnetosphere 7 

particles with ijj -C 9 is still characteristic of pulsar magnetosphere, and in the present paper we concentrate 
on this case of small pitch-angles. 


2.2 Numerical estimates 


As can be seen from the wave equation ( 6 ), the terms containing B can be excluded by introducing the 

substitution E x , y = E x>y e 2c f Bd ~, i.e. these terms act to decrease both field components identically. With 
the common expression for the total intensity, I = E X E X + E y E* (where the asterisk stands for complex 
conjugation), one can easily recognize the coefficient of cyclotron absorption: 


LO 

/x = -QB . 

c 


Using equations (7) and ( 8 ), it can be reduced to the form 


At = 


_, 

mao yto9 2 /2 J 


-f 


( 10 ) 


( 11 ) 


which coincides with th at known for t he coeffic ient of absorption by the system of particles in vacuum (e.g. 
equation (4.13) in iLvubarskii fc Petrovalll998 i. Above it is taken into account that the wave propagation 
is quasi-transverse with respect to the ambient magnetic held, 9 <§C I/ 7 . For the conditions relevant to 
pulsar magnetosphe re, t h e optical deffih to c yclotron absorption, r c = f fJ.dz, is estimated as follows (e.g. 
equation (2.8) in lUviiharskii Petrovalll 99 8?'!: 


r c = 


OAK 2 BH 2 si n 4/5 1 ? 
(P 3 u 971.5 


,3/5 


( 12 ) 


Here re is the plasma multiplicity factor, K 2 = -fp-, B * the magnetic held strength at the surface of the 
neutron star, B* 12 = 1( B ] G , 1 ? the angle between the rotational and magnetic axes of a pulsar, P the pulsar 
period, v the radio frequency, vg = 1 |) j / Hz , 71.5 = 7 s • 0 lle can see that at low enough frequencies the 
absorption depth can be signihcant, especially for the short-period pulsars, P ~ 0.1 s. Below we shall no 
longer be interested in the consequent intensity decrease (equal for both natural waves) and concentrate on 
examining the normalized Stokes parameters with an eye to tracing polarization evolution of the waves. 

In the wave equation ( 6 ), the terms containing G allow for the plasma gyrotropy: if the distribution 
functions of electrons and positrons are identical, G = 0. The problem on the plasma motion in the 
electromagnetic fields of pulsar magnetosphere is very complex and the self-consistent solution has not been 
found yet. Thus, the question on the net current density in the magnetosphere is still open. The simplest 
possibility is that the electrons and positrons differ only in the number densities, with the difference being 
of order of the Goldreich-Julian number density, i.e. A n/n ~ 1/re. Then the relative contribution of the 
plasma gyrotropy to polarization evolution can be estimated as KG/ [{b x + b y ) KA] ~ K~ 1 (9'y) 2 u >'/uoh ~ 
O.lre^ 1 0 ?. 1 7 i 5 io'/ loh- In the resonance region, it can become significant, at least for the rays of a specific 
geometry. However, keeping in mind that true form of G is unknown, we leave out the plasma gyrotropy 
and the resultant cllipticity of the natural waves. Our aim is to study polarization evolution of the linearly 
polarized natural waves and, in particular, to demonstrate how the cllipticity arises and changes purely on 
ac count of wave mode co uplin g and cyclotron absorpt ion. An opposite approach has recently been developed 
bv lLuo fc Melrosel ifeOQdll and lMelrose fc LikJ J2004ali . who investigated the characteristics of the elliptically 
polarized natural waves ignoring the above mentioned effects of polarization evolution. The propagatio n of 
the e lliptical natural waves through the region of cyclotron resonance has been considered in llVTelrose T.i 1 (1 

jhofMbf). 
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8 S. A. Petrova 


In the present paper, we neglect the terms B and G in equation ( 6 ) and concentrate solely on A. Note 
that UM describes wave mode coupling, while QA corresponds to cyclotron absorption. In the resonance 
region, these two contributions are roughly of the same order. 

Taking into account that the magnetic field strength decreases with distance from the neutron star as 
z~ 3 , one can estimate the radius of cyclotron resonance, where u>j9 2 /2 = ujh, for the particles with some 
characteristic Lorentz-factor 7 as 


tl 


= 0.5 5P- 1 Bl / 1 3 2 u- 1/3 1 ~l /3 eZ 2 1 /3 . 


(13) 


Here 77 , = 5 • 10 9 P cm is the light cylinder radius, 6 *_ 1 = 6/ 0.1. The altitude of the mode coupling region, 
z p is another basic quantity of the problem considered. Proceeding from the definition of z p , 

w _ _ , 

7 3 w 2 <9 2 c 

(for more detail see, e.g. lPetrovaIl2003alh one can obtain the following estimate: 


Zp 

tl 


= 0.18P 


-3/2^-3/2 ^1/2 b 1/2 7 -1/2q-1 


7l.5 


*12 "9 


e: 


(14) 


Hence, the altitudes of the regions of pronounced polarization evolution are in the ratio 


Zp 

Z c 


0.33P- 1 ^-l /6 9ll /3 Bl{ 6 2 v- 1 ' 6 K- 1 ' 2 . 


(15) 


One can see that for the parameters relevant to pulsars this ratio can take the values more and less than 
unity. These two possibilities imply qualitatively different scenarios of polarization evolution. In the next 
Section, we dwell on the analytical treatment of the wave equation in the limiting cases z v lz c 1 and 
Zpjz c -C 1 , and find out general features of polarization behaviour in these two scenarios. 


3 ANALYTICAL TREATMENT OF POLARIZATION EVOLUTION 

In terms of the Stokes parameters, 

I — E X E X + E y Ey , 
q — E X E X EyEy , 
u = E X E* + E y E* , 
v = i(E x E* - E y E x ), 

where the asterisk stands for complex conjugation, the wave equations ( 6 ) are written as 
^ = 2 A 2 b x b y u + A 2 (b 2 x - b 2 v )q + A 2 (b 2 x + b 2 y )I , 

-j— = 2Aib x b y v + A 2 (b x — b y )I + A 2 (b x + b y )q , 

^ = 2A 2 b x b y I - Ax(b 2 x - b 2 )q + A 2 {b 2 x + b 2 )u, 

^7 = -2A 1 b x b y q + Ai( 6 ^ - b 2 )u + A 2 {b 2 x + b 2 y )v . (16) 

Here A\ = A 2 = and, in accordance with the above discussion, the terms containing B 

and G are left aside. The last terms in equations (16) can be excluded using the substitution s m = 
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s m exp[— f A 2 (b% + by)dz], m = 1,4, where s m stands for the m-th Stokes parameter. Our primary 
interest is the evolution of the normalized Stokes parameters, so hereafter we consider the quantities s m 
and omit the tildes. 

Introducing the normalized coordinate w = z p /z : one can present equations (16) in the following form: 
d I_ 
dw 
dq_ 
dw 
c\n 

— = wF lV - 2 vF 2 I , 
dw 

dv 

dw 

where 


= —wFiq — 2[iF2U , 
= —wF'i I — 2fj,Fiv , 


= —wF\u + 2fj,Fiq , 


(17) 


F\ = v.p. 


/(7)d7 


(7/7o) 3 [1 - (7/7o) 2 («/« c ) 6 ] 

Fi = fr (. z/zcf 70 / [( z/z c )~ 3 70 ] , 


70 is the characteristic Lorentz-factor of the plasma particles with the distribution function /( 7 ) and z c 
is the characteristic altitude of the resonance region defined by the condition u>h(z c ) = W7o# 2 /2. Due to 
continuity of the plasma flow in the tube of open magnetic lines, the plasma number density n oc B oc z~ 3 . 
The radio waves are emitted in the plane of magnetic lines and hence initially b y = 0. As the waves propagate 
in the rotating magnetosphere, they leave the plane of magnetic lines, so that b y slightly increases along 
the trajectory, b y ~ z/tl- For the sake of simplicity, in the present paper this angle is assumed to remain 
much less than the wavevector tilt in the plane of magnetic lines, b x \ fi = ( b y /b x ) Zp -C 1. 

Strictly speaking, the wave equations ( 6 ) correspond to the frame corotating with the neutron star; for 
the o bserver’s frame the corrections for rotational aberration should be included (cf., e.g. iLvubarskii fz. Petroval 
11999 7 Similarly to b y , these corrections are also ~ zpL and depend on the angle between the rotational and 
magnetic axes of the pulsar, the observational geometry and pulse longitude. However, here we are interested 
only in the principal features of polarization evolution, so we do not specify the geometrical parameters 
and assume that /.i already includes the factor of order unity which allows for rotational aberration. 

It is convenient to assume that the particle Lorentz-factors lie within a finite interval [ 71 , 72 ]- Then 
the resonance region is also finite, with the boundaries given by the conditions uih(z c1 2 ) =w 7 i, 2 # 2 / 2 , and 
the wave trajectory is divided into the three segments lying below, inside and above the resonance region, 
respectively. Let us consider the polarization transfer in each of these regions separately. 


3.1 The region before cyclotron resonance 

Below the region of cyclotron resonance, F 2 = 0, and the set of equations (17) is markedly simplified. For 
the natural waves of pulsar plasma, the initial conditions read: qo/Io = ± 1 , uo = Vo = 0 , where the upper 
and lower signs correspond to the ordinary and extraordinary waves, respectively. It is easy to see that if 
fi <C 1, in the region considered the quantities |u| and |u| remain ~ /u, while |q| « 1. Thus, equation (17) 
can be reduced to the form 
du 

— = wFiv , 
dw 

du 

dw 


= —wF\u + 2/uFj . 


( 18 ) 
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Substituting the second of the above equations into the first one, one can find the solution: 

ui = dz2fi^sRi , v\ = = j = 2/i3^i?i, (19) 


where 

r-, —i f Fiw'dw' I if F-twdLw j-, , / 

Ri = e J / e J Fidw . 

j w 


( 20 ) 


The above solution differs substantially for the cases wi ^$> 1 and wi 1, where wi is the lower boundary 
of the resonance region. It is convenient to introduce the variable x = w/rj = z c /z, which remains of order 
unity, and estimate the integral in equation ( 20 ) at r/ 1 and r/ <C 1 . 

Given that rj^> 1, we have 


mi 



mi 


T 


2fi 

Firfx 2, 


( 21 ) 


which is a common solution in the regime of geometrical optics. In the opposite limit, rj <C 1, the waves 
pass through the region of wave mode coupling. It is reasonable to present the integral (20) as a sum of 
the two integrals over the intervals oo > w > y/rj and y/rj > w > r/ui. The boundary between the intervals, 
u a = 1 / y/rj , w a = y/rj, is chosen so that w is small enough for the wave mode coupling to be efficient only 
within the first interval and at the same time it is far enough from the resonance region. In the first interval, 
the solution has the form 


ui a « ±ny/n(l — w 2 /2) , 

vi a « =F//v^(l + m> 2 /2) ± 2/aui . (22) 


To the first order in rj, for the second interval we have the following solution: 
mi = ±/j,y/n , 

r u i 

mi = i 2/M7 / Fidu ± 2^iy/rj. (23) 

J i/Vv 

Comparing the solutions (21) and (23) obtained for the limiting cases rj 1 and 77 <g 1, respectively, 
one can see that the waves coming to the resonance region can have substantially distinct polarization 
characteristics. 


3.2 The resonance region 


I 11 the resonance region, the set of equations (17) has the form 

d I 2 „ 

- = - V xF 2Q , 

= rfxFiv — 2[ir}F2l , 

da: 

4 ^ = —r) 2 xF\u + 2/j.riFiq , 

da; 

where the terms ~ fj 2 are omitted. For the first two equations, one can find the following solutions: 


q a = Io = exp 


Fha/da/ 


(24) 
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Qe — -le — 6xp 


-q 


F 2 xdx' 


(25) 


where the subscripts ”o” and ”e” correspond to the original ordinary and extraordinary waves, respectively. 
One can see that the two types of waves evolve in a different manner. It should be kept in mind, how¬ 
ever, that actually the ordinary mode intensity does not increase, since it contains an additional factor 
exp[— q 2 ff 1 F 2 x'dx'\ omitted after equation (16). Thus, the ordinary wave is absorbed by the plasma just 
as by the system of particles in vacuum (the absorption coefficient is given by equation (11)), whereas the 
extraordinary wave is absorbed somewhat more efficiently. 

Using equations (25) in the last two equations of the set (24), one can obtain: 


vn = ±—SsRn + vo cos I p I Fix'dx' I — 

q 11 ' 


uo T 


2 q 
qxi 


sin q / Fix'dx 1 


Mn — i——3if?ii± —— exp 


qx 


dhq / F 2 x'dx' 


+wo sin q / Fix'dx' ) + 


Mo =F 


2 /x 
qx i 


cos q / Fix'dx 


where 


Rn = exp 


— iq 2 / Fix'dx' 


exp 


iq 2 f x Fixdx exp ff F 2 xdx 


x' 


—d*'. 


w 0 = Mi(xi), uo = ui(xi). 


In the limit q 1, we are interested only in the evolution of the ordinary wave, since the extraordinary 
one is severely absorbed. The asymptotic treatment of equation (26) yields: 


(«/- 0 n 


2/r 
qx ’ 


(v/I)u = 


2q 

Fiq 3 x 3 


(27) 


This solution still coincides with that known for the regime of geometrical optics in the infinitely strong 
magnetic held (cf. equation (21)). Thus, if in the resonance region the wave propagation obeys geometrical 
optics, i.e. if q 1, the cyclotron absorption does not affect polarization state of the wave and acts only 
to suppress the total intensity. 

In the limit q 1, equation (26) is reduced to the following form: 


Mu = dz/xpiv — 2qq / F 2 dx ', 

J X\ 


rx 

Mil = TMx/tt ± 2qq / F\dx ± 2/j.^/rj. (28) 

ll/vflj 

In contrast to the case of geometrical optics, in case q -C 1 the cyclotron absorption changes the polarization 
state of the waves (cf. mi and mii). This can be understood as follows. The waves are subject to wave mode 
coupling well before the resonance region and become a coherent mixture of the two natural waves. Since 
these constituents are absorbed not identically, the resultant polarization is modified. For the original 
ordinary and extraordinary waves, the amplitude ratio of the entering natural waves is essentially distinct, 
so they suffer different polarization evolution (see equation (28)). 


3.3 The region beyond cyclotron resonance 

As soon as the waves go out of the resonance region, their intensities stop changing, whereas the polarization 
still evolves. The final polarization state of the waves at infinity is as follows: 
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0 f X2 


0 r 2 

Uf = 2/^077111 + uo cos 

r/ / Fixdx 

Jo 

— vo sin 

77 / F\xdx 

Jo 


Vf = —2fLrfRRm + uo sin 

rx 2 

rj 2 / F\xdx 

+ Vo cos 

rx 2 

rf / F\xdx 

where 

Jo 


Jo 


2 r° C X2 o r x ii 

r~. — in * | Fi xdx I iri * I Fi x dx r-, ■, 

Rui = e J / e J Fidx , 

Jo 

u 0 = un(x 2 ), v 0 = vn(x 2 ). 


In case rj 1, equation (29) is reduced to 

_ rx 2 

Uf — uo = ±/Ltv^ — 2 /rr 7 / F 2 dx , 

J X\ 

f° 

Vf = FvV™ ± 2VV / Fidx ± 2ny/rj. 

Ji/Vv 


(30) 


Hence, the wave polarization does not suffer marked changes beyond the resonance region. 

Given that rj 1, in the region under consideration, the wave polarization evolves drastically because 
of wave mode coupling. In equation (29), the addends containing u o and vo compensate for the first terms 
of the expansion of Rui in power of r; -2 , and the main contribution comes from the stationary point of 
Rm. The final polarization takes the form 

( u/I)f = 2 6 / 8 /x?7 -3 ' /4 r(7/8) cos(7r/16), 

(v/I) f = 2 5 / 8 W - 3 / 4 T(7/8) sin(7r/16), (31) 


where T(-) is the gamma-function and it is taken into account that far from the resonance region F\ ss —u 6 , 
independently of the detailed form of the distribution function /(y). 

Equations (30) and (31) describe the final polarization of the waves in cases when the region of wave 
mode coupling lies much lower and much higher than the region of cyclotron resonance, respectively. It 
should be noted that in the weak magnetic field polarization evolution as a result of wave mode coupling 
differs essentially from that in the superstrong magnetic field. First of all, the resultant circular polarization 
has opposite signs (cf. equations (30) and (31)). Besides that, in the weak field, the evolution is less 
pronounced (note the factors jy -3 ' 4 in equation (31)). 


4 NUMERICAL SIMULATION 


Now let us turn to numerical simulation of polarization transfer in pulsar plasma based on equation (16). 
First of all, we introduce a simplified distribution function of the plasma particles in the form of a triangle: 


/( 7 ) 


7 - 7o(l - a) 

-3- 

a 7o 

7o(l + a) - 7 

-^- 3 “ 2 - 

a 7o 


7o (1 - a) ^ 7 ^ 70 , 
7o < 7 < 7o(l + a) ■ 


(32) 


Figure 1 shows the functions F\(z c /z) and F 2 (z c /z), which incorporate equation (32), for the two values of 
the width of /( 7 ): a = 0.1 and a = 0.5. 

According to the approximate solutions (30) and (31), polarization evolution should be most pro¬ 
nounced at rj « 1. In Fig. 2, we present the numerically simulated evolution of the normalized Stokes 
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parameters at 77 = 1 . One can see that the changes are indeed marked and are not identical for the two 
types of waves. 

The final values of the normalized Stokes parameters u/I and v/I as functions of 77 at different 77 are 
plotted in Figs. 3 a and b, respectively. The analytical approximations given by equations ( 30 ) and ( 31 ) are 
shown in asterisks (where appropriate). The polarization parameters orthogonal to those of the original 
ordinary mode are shown by dotted lines for comparison. As can be seen from Fig. 3 , at 77 ~ 1 the modes 
are markedly non-orthogonal. Note also the switch of the sign of v at 77 ~ 1 . 

At large enough 77, the extraordinary mode abruptly turns into the ordinary one. This can be understood 
as follows. Generally speaking, the natural modes are completely independent only infinitely far from the 
region of wave mode coupling, whereas for finite 77 in the resonance region each of the waves contains already 
a small admixture of the other mode. For large enough r /, the extraordinary component is absorbed much 
more efficiently than the ordinary one, so that even in the original extraordinary wave only the ordinary 
component survives. Note, however, that this regime can hardly be observed, because the ordinary wave 
intensity is also strongly suppressed by cyclotron absorption. 


5 DISCUSSION 


As is found above, polarization evolution of radio waves in pulsar magnetosphere can be substantial. The 
wave polarization varies because of wave mode coupling and cyclotron absorption, the change being most 
pronounced if the regions of the coupling and resonance approximately coincide. 

In cases when the resonance region lies well below and well above the region of wave mode coupling 
(77 1 and 77 <g 1 , respectively), the resultant polarization is essentially distinct. First of all, this distinction 

is caused by the different character of the wave mode coupling in the approximations of the weak and strong 
magnetic field. In the former case, polarization evolution is weaker and the sign of circular polarization is 
different. 

For jj > 1 and 77 <C 1 , the consequences of cyclotron absorption are also distinct. Given that 77 1 , the 

waves passing through the resonance region still obey the geometrical optics approximation and propagate 
independently. Then the cyclotron absorption affects only their intensities. In contrast to the cyclotron 
absorption by a system of particles, in the plasma the wave intensities are suppressed not identically: the 
extraordinary wave is absorbed more efficiently, the difference strongly increasing with 77. The numerical 
calculation shows that even at 77 = 1 (/r = 0 . 3 , a = 0 . 1 ) the outgoing intensities differ by a factor of about 
4 . 


In case 77 <g 1 , the cyclotron absorption acts mainly to change the wave polarization. The waves entering 
the resonance region present already a coherent mixture of the two types of natural waves, and since these 
constituents are absorbed not identically, the resultant polarization changes. For the original ordinary and 
extraordinary waves, polarization evolution on account of cyclotron absorption is distinct, so that they 
become non-orthogonal. 

The features of polarization transfer fisted above are believed to have a number of observational con¬ 
sequences. For example, the pulse-to-pulse fluctuations of 77 aro und 77 » 1 can cause the sense reversals 
of the circular polarization at a given p ulse longitude (see also iMelrose fc Fnoll 2 flf) 4 bll . which are really 
observed (e.g. |Karastqrgiou qt alJ | 2 Qq 3 ) ■ It should be kept in mind, however, that this requires too large 
fluctuations in the plasma parameters; besides that, on both sides of 77 = 1, the absolute values of v are 
markedly different, while very often the observed sense reversals do not change the absolute value of the 
circular polarization essentially. 

Generally speaking, pulsar radiation presents a mixture of the two types of waves with nearly orthogonal 
polarizations and comparable randomly varying intensities, so that it is partially depolarized. The cyclotron 
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absorption, especially at large enough 77, introduces a systematic difference in the randomly varying mode 
intensities. Then the ordinary wave should prevail on average, in which case the degree of polarization 
should be markedly enhanced. Large 77 do not seem typical of the pulsars observed, however, some pulsars 
with exceptionally high polarization indeed exist. The Vela presents a well-known example of such a pulsar. 
As a rule, in these cases the same polarization mode dominates throughout the average pulse and in most 
of the single pulses. It is interesting to n ote t h at a few stron gly polarized pulsars have the total-intensity 
profiles classified as partial cones (e.g. Ibvne fc Mauchesteill 988h . which is an additional hint at a significant 
cyclotron absorption in these objects. 

Unfortunately, the actual value of the position angle of linear polarization is not determined in ob¬ 
servations and remains unknown, so that direct observations tell nothing as to whether the ordinary or 
extraordinary mode dominates in the strongly polarized pulsars. Even in case of the Vela, where the X- 
ray observations of the wind nebula give some new insights into the pulsar geometry, this question has 
no firm answer. As for the theoretical predictions, cyclotron absorption strongly favours the dominance of 
the ordinary waves. Let us note in passing that in case of resonant absorption by the particles with large 
pitch-angles , which is a very probable regime for the Vela, the ordinary waves should also prevail llPetroval 

boodhoo.m 

As can be seen from equation ( 15 ), 77 is a very weak function of radio frequency. Obviously, just for this 
reason the majority of pulsars are detectable over a wide frequency range. Let us speculate, however, that if 
at the intermediate frequencies 77 ~ 1, at low frequencies 77 is larger and the degree of polarizati o n can sub¬ 
stant ially increase, which is in line with the observational results below 100 MHz liSulevmanova V Pugachevl 
hoodi . At high frequencies 77 < 1 and the polarization modes can be markedly non-orthogonal, in which case 
the degree of circular polarization can be enhanced. The increase of circu lar polarization at h igh frequencies 
is indeed found in the average profiles of a number of pulsars (e.g. Ivon Hoensbroech fc Lesc 3 S). 

A possible non-orthogonality of polarization modes on account of cyclotron absorption has a number of 
important implications for the statistical model of the individual-pulse polarization of pulsars. (This issue 
will be considered in detail in the forthcoming naner IPetrovalKoOft ). Given that the observed superposed 
orthogonal polarization modes are associated with the natural waves of pulsar plasma, they should be 
orthogonal by definition. The process of wave mode coupling also preserves the orthogonality of the waves. 
To the best of our knowledge, cyclotron absorption is the only way for the superposed modes to become 
non-orthogonal. Note, however, that the non-orthogonal states of the observed sum of the modes in different 
single pulses can be attributed to the time-dependent action of propagation effects. 


6 CONCLUSIONS 

We have studied polarization transfer in the hot magnetized plasma of pulsars. In the present consideration, 
we have assumed the non-gyrotropic plasma, with the identical distributions of the electrons and positrons, 
and the small pitch-angles of the particles. Proceeding from the Maxwell’s equations, we have derived the 
set of equations describing the evolution of the Stokes parameters of the original linearly polarized natural 
waves. These equations have been solved analytically and the results have been confirmed by numerical 
calculations. 

The polarization evolution of the waves has been found significant. The polarization characteristics 
change on account of the wave mode coupling and cyclotron absorption. In cases when the region of coupling 
lies well above and well below the resonance region, the resultant polarization is qualitatively distinct. If the 
waves pass through the region of wave mode coupling first, they acquire the elliptical polarizations purely 
orthogonal at the Poincare sphere. Further on, in the resonance region, they become non-orthogonal. If the 
waves enter the resonance region before coupling, cyclotron absorption does not affect their polarization 
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states and suppresses the total intensities only, the extraordinary wave being absorbed somewhat more 
efficiently. Further on the waves suffer mode coupling in the limit of weak magnetic field, which differs from 
that in the strong field. Firstly, the total change of polarization parameters is much less. Besides that, the 
resultant circular polarization has the opposite sense. 

The observational consequences of polarization transfer in pulsar magnetosphere can be summarized 
as follows. Because of cyclotron absorption, at large enough r/ one mode can markedly dominate another 
one, in which case one can expect strongly polarized profiles, with the same mode dominating throughout 
the average pulse and in most of the individual pulses. This is indeed characteristic of several pulsars. In 
some cases, the total-intensity profiles of the strongly polarized pulsars are classified as partial cones, with 
one of the conal components being absent. This seems to be an additional argument in favour of strong 
cyclotron absorption in these pulsars. At present it is not known exactly what of the natural waves actually 
dominates in the strongly polarized pulsars. Our consideration favours the dominance of the ordinary mode. 

If one consider i) as a function of radio frequency, with 77 « 1 at the intermediate frequencies, in the 
low-frequency range r] ^ 1 and one can expect strong polarization of pulsar profiles, whereas at high enough 
frequencies there may be an increase of the degree of circular polarization because of non-orthogonality of 
the modes. The non-orthogonality of the outgoing waves is of a crucial importance for the statistical model 
of the individual-pulse polarization. This will be studied in detail in the forthcoming paper. 
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Figure 2. Evolution of the normalized Stokes parameters of the original ordinary (a-c) and extraordinary (d-f) 
waves along the trajectory; r) = 1, /x = 0.3, a = 0.1. 
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Figure 3. The final normalized Stokes parameters as functions of 77 for different /a. The analytical approximations 
given by equations (30) and (31) are shown in asterisks. The solid and dashed lines correspond to the original 
ordinary and extraordinary waves, respectively. The dotted lines show the parameters orthogonal to those of the 
ordinary mode. 
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APPENDIX A: CONDUCTIVITY OF THE HOT MAGNETIZED PLASMA 


Our aim here is to find the conductivity tensor of the hot magnetized plasma in the coordinate system 
with the z-axis along the wavevector. In the commonly used system with the z-axis along the magnetic field 
direction, it has the following form (e.g. 




Here the brackets < ... > stand for the integral operator J ...p±dp±dp z ; J n and J' n are the Bessel function 
of the n-th order and its derivative, £ is their argument, 

£ = k±v±y/ui H ; 

Cn = (w - k z v z - nwH/y )- 1 ; 



v± dp± u> l dp z v± dp± J 


v± dp± to 


1 dF + K 


1 dF + k z f dF v z dF 


_ dF nuj H f dF v z OF \ 

11 dp z yui y<9p z v± dp a. ) ’ 

F is the particle distribution function with the normalization 

J F(p± ,p z )p±dp±dp z = N 0 , 

and the subscripts ”z” and ”_L” refer to the components parallel and perpendicular to the magnetic field, 
respectively. 

In terms of the quantities in the system with Oz||k, v z = U||, k z = kb z , k± = kb±, b± = yrni + bf. 
Then the argument of the Bessel function is written as 

kb±v j _7 uj' ip 

ojh u>h b _L 

where u>' = uiy(b^. + by) is the frequency in the rest frame of the particle, ip the particle pitch-angle and it is 
taken into account that kv/oo sa 1, i.e. v ~ c and the refractive index is approximately unity. Hereafter it is 
assumed that the particle pitch-angle is small, ip ^ b ±, in which case £ < 1 up to the cyclotron resonance, 
a/ = u>h- As can be seen from the final equations ( 6 ), beyond the cyclotron resonance, at u>' > u>h, 
the evolution of the electric field amplitudes rapidly ceases because of strong decrease of u>h, u>h oc r~ 3 . 
Therefore in the case of interest £ <C 1 and one can use the approximation of the Bessel function: 



Then the form of the conductivity tensor is substantially simplified: 




© 0000 RAS, MNRAS 000, 000-000 













22 S. A. Petrova 


. 2 /7^6±/?(l - f3b z )\ 
a ° 13 ~ CT ° 31 _ 16 \ m (, - w ' 2 ) / ’ 


(A2) 


^023 ^032 ^ 


LOHb±f3 
m (uijj — u}' 2 ) 


a 033 — 


'yu!b 2 L (3 2 


\ (uijj — u>' 2 ) m'yui' 2 / 


where the brackets stand for the integral operator J ...F{p\\ )dpy. 
In the system with 0«||k, the conductivity tensor is given by 

<7 = RaoR 1 , 


(A3) 


where R and R 1 are the matrices of the direct and inverse transformation between the coordinate systems 
considered: 


R = 


R- 1 = 



( b ± 


by 
b± 

V b x 


b± 

by 


(A4) 


(AS) 


As is discussed after equation (5), the waves are almost transverse and we take E z = 0. Then we should 
find only 


jx — (JxxEx T (JxyEy and jy — CfyxEx A 0~yy Ely . 

Substituting equations (A4) and (A5) in equation (A3) and taking into account that ao 1[L 
—< 7 o 21 , oo 13 = ( 7 o 31 and <ro 23 = —co 32 , one obtains: 

_ t 2\ . ,2 2& 2 &, 

O'xx 0 X J O'On 4~ O x (J O33 ^ ^0i3 ? 

i2\ . ,2 “ibybz 

®yy Oy) cron by(T o 33 ^ ^013 ? 

crxy — bxby ( + 0*033 2(^2 /b_L )0*013) + b z cro 12 bj_(J02 3 > 

0*2/a: b x by ( cr 0ll + 0033 2(6 2 /6 _l)(To 13 ) ^tz0t)i2 ^_i_ 0 "O 23 5 

Using equation (A2) in equation (A7) yields finally: 

G XX b x Oj2 1 oyy CL\by Q>2 1 

&xy — CL\b x by + Q , (Tycc = CLl_b x by Q , 


with 


ai 


• 2 2 
£e ojuj h 


^(ojI-u' 2 ) 


CL2 


-ie 2 7 aj(l - f)b z ) 2 ) \ 
m (ai|f — a/ 2 ) / 


(A6) 


^022 1 ^0l2 


(A7) 


(AS) 
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g = 


e 2 uj H (b z — 0) \ 
m (uijj — ui' 2 ) J 


Throughout this section, we have considered the one-component plasma for simplicity. The results can be 
extended to the case of pair plasma by introducing the summation over the particle species. 
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